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1. Introduction 

Since there are an equal number of left-handed and right-handed particles in the lattice fermion 
propagator for each set of quantum numbers, there is no chiral anomaly on the lattice for the 
naive fermions. Canonical local formulations, such as the Wilson fermions, the Kogut-Susskind 
fermions, the twisted mass fermions or the Creutz-Borici fermions break the flavour singlet axial 
U(l) and some or all of the flavour non-singlet axial symmetries. The overlap fermions, on the 
other hand, preserve all chiral symmetries and have an index theorem as well. 

Interestingly, chiral anomaly plays a crucial role in deciding the structure of the QCD phase 
diagram for our world with two light quarks and one moderately heavy quark. A fundamental 
aspect of the QCD phase diagram may be the existence of a critical point in the T-}Xb plane, where 
pL B denotes the baryonic chemical potential. The critical point is expected on the basis of chiral 
symmetries and model considerations. Since we employ lattice techniques to investigate the phase 
diagram, the presence on those chiral symmetries on the lattice is desired, if not required, in the 
entire T-pLs plane. The popular staggered fermions have a remnant chiral symmetry on the lattice 
but a not so well-defined flavour number; two flavours are simulated using a square root of the 
fermion determinant for them. Clearly, use of overlap fermions with exact chiral symmetry on 
lattice is desirable. 

Introduction of chemical potential for the overlap fermions turns out to be nontrivial. Ideally, 
one should construct the conserved charge, N, as a first step and then simply add /J.N to the action. 
For the local fermions, this leads to aT 1 divergences in the continuum limit in thermodynamic 
quantities such as the energy density or the number density. One uses the freedom of adding 
irrelevant terms to cure them by multiplying gauge links in positive/negative time direction by 
exp(fl , /i) and exp(— a\i) respectively or more general functions which satisfy certain conditions 
[j|]. Note that this modification, indeed even the simple addition of \iN does not change the chiral 
invariance of the lattice action in anyway. The non-local nature of the overlap fermions makes this 
construction difficult, possibly even non-unique [||]. It was proposed by Bloch and Wettig [^] that 
the same prescription of multiplying the gauge links, as above, may be used. This required them to 
extend the definition of the sign function employed in the overlap operator since the argument of 
the sign function, namely the Wilson-Dirac matrix, was no more 75-hermitian. 

Banerjee, Gavai and Sharma [Q] showed analytically that this prescription does not have any 
aT 1 divergences in the continuum limit but unfortunately the resultant overlap fermion action has 
no chiral invariance on the lattice for nonzero \l. Of course, it is restored in the continuum limit 
but so it is for the simpler local Wilson fermions: in both cases the bare parameters of the action 
such as the quark mass or the chemical potential will get renormalized and depend on the coupling 
g. Alternatively, one can modify the chiral transformation, again by terms which vanish in the 
continuum limit, to restore the chiral invariance. This explains the results obtained by Bloch and 
Wettig that one even has an index theorem for the /I -dependent overlap Dirac operator: the anomaly 
equation itself is jj, -dependent Our work presented in this talk was sparked by the curiosity 
about the fate of the anomaly relation in the continuum on the introduction of nonzero chemical 
potential. 

Before presenting our work let us remark, however, that the change in the chiral transformation 
on the lattice implied by the Bloch-Wettig proposal is undesirable since it makes the generators of 
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the transformation to be non-hermitian. Worse still, the symmetry group itself changes for each and 
every change of making this avenue not very useful for the original problem of investigation of 
the QCD phase diagram in the T-/J.B plane. In order to appreciate the loss, recall that the symmetry 
group remains the same at each T for pi = 0. Therefore, a study of the temperature dependence of 
the order parameter (yy), and any abrupt change in it at the transition temperature, can reveal the 
chiral symmetry restoring transition, if any, and thereby the change in the nature of the vacuum. 
Were the symmetry group to change as a function of T, it would be impossible to attribute a vanish- 
ing chiral condensate uniquely as a change in the vacuum or a phase change, thereby eliminating 
the role of the chiral condensate as the order parameter. Thus one can indeed use the freedom of 
adding terms irrelevant in the continuum limit to modify the chiral transformation on the lattice but 
such a change should not rob one off the usual order parameter for studying the chiral transition. 



2. Anomaly at T = and /i ^ in continuum 

Classically y' = expO'a/s) y and y' = y exp (10175) is a symmetry for the QCD action for 
massless quarks , leading to the current conservation equation d^J^ = 0. At finite temperature 
and/or density, this classical symmetry remains intact. Quantum loop effects can, and do, cause 
corrections already at T = jl = 0, leading to the anomaly equation. One needs to compute (d^J^) 
to check this. The famous calculation of the Adler-Bell-Jackiw(ABJ) triangle diagram for the 
U(l) case []|] demonstrated perturbatively that the axial E7(l) is broken by quantum effects while 
Fujikawa ^ provided a new insight by showing that the anomaly arises due to the change of the 
fermion measure under the corresponding transformation of the fermion fields in the path integral 
method. We have used both the methods to check the fate of the anomaly equation at finite density, 
i.e, on introduction of a nonzero chemical potential jj. at zero temperature. 



2.1 Perturbative calculation 

The lowest order diagram is the canonincal ABJ axial vector-vector-vector(AVV) triangle di- 
agram. It is well-known that the higher order diagrams do not contribute to the anomaly equation 
at zero density, neither do other diagrams like the square and pentagon diagrams. We too therefore 
compute only AW triangle diagram at finite density. Our notation for the QCD Lagrangian in the 
Euclidean space with the finite number density term is the same as in [^]. In order to maintain 
consistency with the lattice literature, we have, however, chosen the Dirac gamma matrices to be 
Hermitian, leading to the action, 
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(2.1) 



where = J v (d v — igA"T a ) with T a being the generators of the SU(3) gauge group. The 75 = 
7i72/374 is also hermitian in our case. Using the canonical Euclidean space Feynman rules, the 
amplitude of these AW triangle diagrams can be computed. Denoting by A^ pa (kife) the total 
amplitude and contracting it with qx, the (d^jn) for fi 7^ can be obtained from the triangle 
diagrams for \i ^ 0: 
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(2.2) 



An inspection of the eq.(|2.2[) reveals quadratically divergent integrals which differ from the 
usual case by merely having /I / 0. The computation therefore can be done the same way by 
writing § q k A^ a = (-/) tr [T a T b ]g 2 J £fr [f{p-k u k 2 ) ~ fipM + /(p-*2,*l) ~ f(p,h)} 
for a suitably chosen function / and introducing a gauge-invariant cut-off. The nonzero /I appears 
in / in the denominator as [(p4 — in) 2 + p 2 ] with linear in /I and a /I -independent term in the 
numerator; the \x 2 term vanishes being proportional to Tr [?????] ~ e 4a4p . The final result 
turns H out to be the same anomaly relation as for /I = 0, since the jx -dependent terms appear with 
A -1 , and vanish as the cut-off A — > oo. The same computation is easily generalized |§] to nonzero 
temperatures. 

2.2 Non-perturbative calculation 

Fujikawa [||] taught us how to compute the chiral anomaly non-perturbatively using the path 
integral formalism. Under the chiral transformation of the fermion fields, the measure changes as 

@\jjr @y = @\jr ^yDet[ ^/ | = ®\jf @y exp(-2/a J d 4 x Try 5 ) . (2.3) 

The trace can be computed using the the eigenvectors of the operator ip. It is an anti-Hermitian 
operator for /I = with purely imaginary eigenvalues and the corresponding eigenvectors form a 
complete orthonormal basis. Using the fact {75,^} = 0, it is easy to show[|8|] that Tr 75 = for the 
space of eigenvectors with nonzero eigenvalues, leading to the usual Tr y 5 = n + — n_ relation. 

ft(n) still anti-commutes with 75 but has both an anti-Hermitian and a Hermitian term. Re- 
markably, it turns out [||] still to be diagonalizable with right and left eigenvectors which to- 
gether form a complete set using which essentially the same argument as above leads to the same 
Tr 75 = n + — ri- for /I 7^ as well. Note that the zero modes are still defined with respect to 
ip{\l =0). This should be contrasted with the index theorem for the overlap operator proposed 
by Bloch-Wettig where the zero modes are jl-dependent [^]. This is perhaps natural since the 
corresponding chiral transformation is jj,-dependent, as discussed above. 

Curiously, a "gauge-like" symmetry, defined by a non-unitary transformation of the fermion 
fields, given by y (x,t) = exp (/it) 1/a(x, t), ^'(x, t) = V^(x,t) exp(— /it) , eliminates all the /l- 
dependent terms of the QCD action. This transformation commutes with the chiral transformations, 
explaining thus the preservation of the same anomaly relation for /I 7^ 0. If one were to insist on 
this symmetry on the lattice, then addition of a simple /J.N is easily seen to be forbidden for any 
local fermion action. Moreover, it can be easily implemented for any such local action, and leads 
naturally to all the known proposals ^ of introduction of the chemical potential on the lattice so 
far. 

While the symmetry clearly arises in the continuum due to the local nature of the fermion 
action, one may choose to demand it for even non-local lattice actions. A natural proposal then 
for the overlap operator for nonzero jj. is D ov (/i) = exp (—/it) D ov (/i = 0) exp (/it). While this 
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can be shown to i) satisfy the Ginsparg- Wilson relation [^] and ii) lead divergence-free ideal gas 
results in the continuum limit, it clearly does not commute with the nonlocal chiral transformations 
unless they are made /i -dependent too. Its main advantage though could still be that the usual sign 
function is used in this overlap operator. 

3. Two simple ideas for the lattice QCD at finite density 
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Figure 1: The energy density (left panel) in units of T 4 and the quark number susceptibility (right panel), 
normalized to the continuum ideal gas value, as a function of 1 /Nj for M values as indicated for lattices 
with aspect ratio four. 

Spurred by the results of the previous sections, we propose two simple ideas for simulations of 
lattice QCD at finite density. Let us state them at the outset: 1) If chiral invariance and the anomaly 
relation are to be only restored in the continuum limit, one may introduce the chemical potential on 
the lattice in a simpler way and 2) known techniques in the literature to eliminate the free theory 
related divergences may simplify the computations. 

From the well-known relation between the domain wall fermions and the overlap fermions 



[ |10| , [lip , we know that only fermions confined to the domain walls are physical. Introducing a 
chemical potential only to count them, one has 



£U£)xy = (D. 



ov )xy 



afX 
2a 4 M 



(1 - y 4 )U 4 (y) 8 xy _ 4 + (74 + l)!^*)^ 



(3-1) 



The chief advantage of eq.(gjj) is that D ov is defined by the usual sign-function for a Hermi- 
tian matrix, making the computations simpler. Of course, it too breaks chiral invariance at the 
same order as the Bloch-Wettig proposal. One has, however, to expect a~ 2 -divergences as a — > 0. 
Following the same prescription which is used for the pressure computation (which diverges at 
zero temperature as A 4 ), we use physical quantities computed on large Nj and the same lattice 
spacing a for subtraction of these divergences. We tested this by considering two physical quan- 
tities for the ideal gas, namely, the difference in the energy density due to nonzero jit, defined by 
Ae(n,T) = e(n,T)-e(0,T) and the quark number susceptibility, % = T /V d 2 ln3f/dn 2 . 

Figure |I] shows the results for these quantities, where we used n/T = 0.5 for the Ae and pL = 
for the x to compare with our earlier [ffl] results using the exponential prescription of Bloch-Wettig. 
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Figure 2: The variation of fourth order susceptibility, normalized by its corresponding continuum value, as 
a function of 1 /Nj for £ = 4, M = 1 .5 for A) ft /s and B)/i/M way of incorporating the chemical potential. 



The zero temperature values were computed numerically on a lattice with a very large temporal 
extent Nt and fixed an, such that T = 1/ (Nja^ — > 0. The Matsubara frequencies then become 
continuous and hence could be integrated upon numerically. The M-values are as indicated along 
the symbol used. From a comparison of the plots with the corresponding ones [Q] for the Bloch- 
Wettig case, we find that i) the zero temperature subtraction procedure does indeed eliminate the 
divergences, ii) there are no oscillations for odd-even values of Nj, iii)the M-dependence is much 
less pronounced, and iv) the scaling towards the continuum value is linear with the possibility of 
an easier extrapolation. 

A further test of the subtraction procedure is that no additional divergences should be seen in 
other quantities as well, e.g., the fourth order susceptibility. Figure ||] demonstrates this to be the 
case. The two curves shown are different possible ways of normalizing the chemical potential term, 
as indicated. Further improvements in the form to achieve a faster convergence to the continuum 
needs to be explored. 

Taking a cue from eq.(|3.1[), one can introduce /I in general by 

S F = Y*{x)M(ji;x,y)V(y) = ZV(x)D(x,y)V(y) + fla^N^y) , (3.2) 

x,y x,y x,y 

where D can be the the staggered, the Wilson-Dirac or any other suitable fermion operator and 
N(x,y) is the merely the corresponding point-split and gauge invariant number density. Clearly, 
as above, one ought to be able to get rid of the annoying divergences by subtracting the same 
physical quantity computed on large Nj and the same lattice spacing a. Usually, one would in- 
troduce exp(±a/i) factors in M and not add the N term. Denoting by superscripted primes var- 
ious derivatives with respect to a\l, one sees that the usual case has M' = M'" ... = Y<x,yN(x,y) 
and M" = M"" = M""" ... ^ 0, i.e, all derivative terms contribute making the successive Tay- 
lor coefficients have more and more terms. In contrast, eq.(3.2) leads to M' = Y.xyN(x,y), and 
M" = M'" = M"" ... = 0, i.e., a lot fewer terms in the Taylor coefficients. This could potentially 
help in obtaining higher order coefficients than computed so far. For example, the 4th (8th) order 
susceptibility, has 4th (8th) derivative of M, which has only one (one) term for eq. (3.2) in contrast 
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to five (eighteen!) in the usual case [12]. Whether the subsequent reduction in the computations of 
M _1 actually helps in speeding these computations is currently being investigated by us. 

4. Summary 

We showed, both perturbatively and non-perturbatively, that the introduction of nonzero chem- 
ical potential, pi, leaves the anomaly unaffected. The zero modes of the Dirac operator for /i = 
govern it. Nonzero pt simply scales the eigenvectors, makes the right and left eigenvectors distinct 
but together they still satisfy a completeness relation. We pointed out further that the reason for 
this can be traced to a "gauge-like" symmetry in the continuum using which /I -dependent terms 
can be "gauged" away. All the currently known prescriptions of introducing pt on lattice can be 
understood in terms of a similar symmetry on the lattice, and ought to protect the anomaly if the 
fermion regularization itself has it. While such a symmetry can be shown to exist for local actions, 
it could be simply adopted for the nonlocal cases such as the overlap fermions. The resultant effec- 
tive overlap operator for nonzero pL also satisfies the Ginsparg- Wilson relation but unfortunately is 
not invariant under the chiral transformation. 

We also showed that investigating the overlap fermions at finite density by simply adding 
the /I -term linearly is feasible. It too has similar chiral symmetry breaking as the Bloch-Wettig 
proposal but the corresponding inverse propagator is simpler and better-defined. Extending the 
idea of addition of pi linearly to the usual staggered fermion case was shown to open a possible 
avenue for computations of higher order coefficients. 
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